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Abstract 

e + e~ annihilation data recorded by the JADE detector at PETRA were used to 
measure the C-parameter for the first time at yfs = 35 and 44 GeV. The distributions 
were compared to a resummed QCD calculation. In addition, we applied extended 
resummed calculations to the total and wide jet broadening variables, Bt and Byy. 
We combined the results on a s with those of our previous study of differential 2- 
jet rate, thrust, and heavy jet mass, obtaining o s (35 GeV) = 0.1448 j^iooTO 

and 

a s (44 GeV) = 0.13921qqq74. Moreover power corrections to the mean values of the 
observables mentioned above were investigated considering the Milan factor and the 
improved prediction for the jet broadening observables. Our study, which considered 
e + e~ data of five event shape observables between yfs = 14 and 183 GeV, yielded 
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1 Introduction 



In a recent publication 0] we have presented a study on event shape variables and deter- 
minations of a s using data of e + e~ annihilations at y/s = 22 to 44 GeV recorded with 
the JADE detector [[| at the PETRA collider. This study provided valuable information 
which was not available before from e + e~-annihilations in the PETRA energy range. The 
results on a s , obtained in a manner similar to the determinations at LEP, demonstrated 
that the energy dependence of a s (Q) is in good agreement with the prediction of Quantum 
Chromodynamics (QCD). Evolved to the Z° mass scale, the results agree well with those 
obtained at LEP, and are of similar precision. In addition, power corrections, applied to 
analytic QCD calculations of the mean values of event shape distributions, were found to 
qualitatively and quantitatively describe the effects of hadronisation. Thus QCD could be 
tested without the need of phenomenological hadronisation models. 

Meanwhile the perturbative calculations for the jet broadening variables were improved 
by including a proper treatment of the quark recoil [|J. Furthermore, a resummation of 
leading and next-to-leading logarithm terms to all orders of a s (NLLA) became available 
for the C-parameter ||. Besides these advances in the perturbative description of event 
shape observables progress was made in the understanding of non-perturbative power 
corrections to the event shape observables and their mean values. In particular two- 
loop calculations of such corrections were performed which modify the one-loop result of 
the power correction to the event shapes by a factor (Milan factor) ||. Revisiting the 
interdependence between perturbative and non-perturbative contributions for the case of 
jet broadening observables yielded in Reference a further modification of the results 
given in [[|. 

In this paper we complement our previous publication by a new 0(a^)+NLLA de- 
termination of a s from C-parameter and update our a s determination from jet broadening 
at \fs = 35 and 44 GeV in the Sections |3] and f|. We also applied power corrections to 
the C-parameter distributions and re-investigated those for the mean values of the thrust, 
heavy jet mass and both jet broadening observables in Section |5|. We start in Section 
with a brief summary of the data samples used and draw conclusions from our results in 
Section || 

2 JADE data 

The JADE detector was operated from 1979 until 1986 at the PETRA electron-positron 
collider at centre-of-mass energies of yfs = 12 to 46.7 GeV. A detailed description of the 
JADE detector can be found in 0], |3[]. The components mainly used for this analysis were 
the central jet chamber to measure charged particle tracks and the lead glass calorimeter 
to measure energy depositions of electromagnetic showers, which both covered almost the 
whole solid angle of Air. 

For the present study we consider data collected by the JADE detector in 1984 to 1986 
at centre-of-mass energies between 39.5-46.7 GeV and at 35 GeV. Multihadronic events 
were selected by the standard JADE selection cuts || which have been summarised in 
detail in our previous publication 0. Background from two-photon processes and r-pair 
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events and from events with hard initial state photon radiation were mainly suppressed 
by posing constraints on the visible energy, the total momentum and the charged particle 
multiplicity of the event, and by cuts in the polar angle of the thrust axis. The estimated 
purity of the resulting data samples is about 98% |§. The numbers of events retained 
after these cuts are 6158 around yfs = 44 GeV and 20 926 at yfs = 35 GeV. 

Corresponding original Monte Carlo detector simulation data for 35 and 44 GeV were 
based on the QCD parton shower plus string fragmentation implemented by the JETSET 



6.3 event generator ||10|| . Both samples included a detailed simulation of the acceptance 
and resolution of the JADE detector. We used the standard set of parameters for the 
event generation as described in f2|. As pointed out there, the agreement of the simulated 
data with the measurements is good. 



3 Measurement of event shapes 

From the data passing the selection criteria as described in the previous section, the 
distributions of the C-parameter and the jet broadening event shape variables Bt and By/ 
were determined using the 3-momenta pi of the reconstructed particles. The C-parameter 
is defined in Reference iTTll as 



C = 3(AiA 2 + A 2 A3 + A 3 Ai) 
where A 7 , 7 = 1,2,3, are the eigenvalues of the momentum tensor 







The definition of the total and wide jet broadening measures can be found in Reference W2 



Limits of the acceptance and resolution of the detector and effects due to initial state 
photon radiation were corrected by a bin-by-bin correction procedure. The correction 
factors were defined by the ratio of the distribution calculated from events generated by 
JETSET 6.3 at hadron level over the same distribution at detector level. The hadron level 
distributions were obtained from JETSET 6.3 generator runs without detector simulation 
and without initial state radiation, using all particles with lifetimes r > 3 • 10~ 10 s. Events 
at detector level contained initial state photon radiation and a detailed simulation of the 
detector response, and were processed in the same way as the data. 

The data distributions were further corrected for hadronisation effects by applying 
bin-by-bin correction factors derived from the distributions at parton level and at hadron 
level, where the parton level is given by the partons before hadronisation. In the upper 
part of Fig. |], the hadronisation correction for the distributions of the C-parameter at 
y/s = 35 and 44 GeV are shown. The correction factors for all event shapes are typically 
of the order 10 to 30% increasing towards the 2-jet region. 

Systematic uncertainties of the corrected data distributions were investigated by mod- 
ifying details of the event selection and of the correction procedure. For each variation 
the whole analysis was repeated and any deviation from the main result was considered 
as a systematic error. In general, the maximum deviation from the main result for each 
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kind of variation was regarded as a symmetric systematic uncertainty. The main result 
was obtained using the default selection and correction procedure as described above. 

To estimate experimental uncertainties, we considered either tracks or clusters for 
the measurement of the event shape distributions. Furthermore, the selection cuts for 
multihadronic events were either loosened or tightened as described in 0, in order to 
check the influence of background events on the measurement. 

The impact of the hadronisation model of the JETSET 6.3 generator was studied by 
varying the values of several significant model parameters around their tuned values from 
Reference (13| used for our main result. The variations amount to the percentage of the one 
standard deviations found in Reference [14] from a parameter tuning of JETSET. In detail, 
effects due to parton shower, hadronisation parameters, and quark masses were considered, 
following the procedure in 0. These uncertainties of the hadronisation corrections are 
shown for the C-parameter distributions at y^i = 35 and 44 GeV as shaded band around 
the correction factors in the upper part of Figure |l[ 



4 Determination of a* at \fs = 35 and 44 GeV 



The data values for the event shape distributions and the corresponding mean value of 
the C-parameter corrected to hadron level are listed in Table |l], where also the statistical 
errors and experimental systematic uncertainties are given. Our measurements of the jet 
broadening event shape distributions can be found in Reference 0. 

After correcting for hadronisation, the event shape distributions were compared directly 
with the analytic QCD calculations. We determined the strong coupling constant a s by 
X 2 fits to the event shape distributions of C, B T , and B w . The fits to the C-parameter 
distributions considered the resummation results obtained in 0. For the fits to the jet 
broadening measures we used the improved calculation of Reference [|J. For the sake of 



direct comparison with other published results we chose the ln(i?)-matching scheme |i~5 
to merge the 0(a^)[|l6] with the NLLA calculations. The renormalisation scale factor, 
x^ = fi/ t/s was set to — 1 for the main result. 

The fit ranges for each observable were chosen to be the largest range for which the 
hadronisation uncertainties remained less than 15%, for which the contribution of an 
extreme bin did not dominate the x 2 /d.o.f., and for which fits were stable. The remaining 
changes in a s when modifying the fit range by one bin on either side were taken as 
systematic uncertainties. Only statistical errors were considered in the fit. 

In order to investigate the importance of higher order terms in the theory, we also 
changed the renormalisation scale factor in the range of x^ = 0.5 to 2.0. The associated 
renormalisation scale uncertainties are larger than those from the detector correction and 
the hadronisation model dependence. 

Our results for a s from Bt, B\y and the C-parameter at = 35 and 44 GeV are given 
in Table [| where also the statistical and systematic uncertainties of the measurements 
are given. It should be pointed out that the improved perturbative calculation for the jet 
broadening resulted in an increased a s value, not affecting the systematic uncertainties 
but improving the x 2 /d.o.f. 

In Figure ^] the a s values from this study are shown with those obtained from the 



4 



distributions of thrust, heavy jet mass and the differential 2-jet rate using the Durham 
scheme as presented in our previous publication 0. Replacing the ot s values obtained 
from the jet broadening observables, a single a s value was obtained from the individual 
determinations from the six event shape observables following the procedure described in 



References |I7[]. This procedure accounts for correlations of the systematic uncertainties. 
At each energy, a weighted average of the six a s values was calculated with the reciprocal 
of the respective squared total error used as a weight. For each of the systematic checks, 
the mean of the a s values from all considered observables was determined. Any deviation 
of this mean from the weighted average of the main result was taken as a systematic 
uncertainty. 

The final results for a s are 

a s (35 GeV) = 0.1448 ± 0.0010(stat.) ^^(syst.) 
a s (44GeV) = 0.1392 ± 0.0017(stat.) +°;°^(syst.) , 



where the systematic errors at 35 and 44 GeV are the quadratic sums of the experimental 

h0.0070 
-0.00355 



uncertainties (±0.0017, ±0.0032), the effects due to the Monte Carlo modelling '+ - 0070 



-aoo2?) an d the contributions due to variation of the renormalisation scale (io!oo57> ^aooss) 
The modelling uncertainties due to quark mass effects contribute significantly to the total 
error. 



5 Power corrections to the mean values of event shapes 

The strong coupling constant a s can also be extracted from the energy dependence of 
the mean values of event shape distributions. Non-perturbative effects are inherent in 
the event shape observables due to the contributions of very low energetic gluons and of 
hadronisation. For an observable T the effect of the gluons can be represented by additive 
power-suppressed corrections ((jF pow ) oc 1/y/s) to the perturbative predictions of the mean 
values of the event shape distributions ((jFP ert )). In the calculations of Reference || [1 
which we used in this analysis a non-perturbative parameter 

1 rvi 

a (/x/) = — / dk a s (k) 



Hi Jo 

was introduced to replace the divergent portion of the perturbative expression for a s (y/s) 
below an infrared matching scale /ij. The power corrections to event shape observables 
have been calculated in References up to two-loops. It was found that the two-loop 



result is identical to the one- loop calculation of Reference |L8| up to a factor M. = 1.79 ± 
0.36, which is known as the Milan factor, and a factor 2/rr, which is due to defining the 
non-perturbative parameter a as above. 

The power corrections for the mean values of thrust (T), heavy jet mass (Mfj/s) and 
C-parameter is 

<^ pow -> = a F ■ f-— ■ • (aoW - a s (^) + 0{a 2 s )) . 

3rr 7T Vv 5 / 
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It was recently found in Reference [0] that an additional factor has to be applied in the 
case of jet broadening observables. In leading order the factors are given by 



7T 7T 

^=^= - 0.822 , — = - 0.343 



4^(0.472^/3 4 v /2a s (0.472 v /i)/3 

for total (B T ) and wide (B w ) jet broadening, respectively. The coefficients which 
depend on the observable T are calculable. They assume the values —1, 1/2, 1/2, 1/4, 37r/2 
for thrust, heavy jet mass, total and wide jet broadening, and C-parameter, respectively || 

The power correction to the mean values of an event shape T is an additive term to 
the perturbative prediction, (J 7 ) = (jF pcrt ) + (jF pow -). We determined a s (M z o) and the 
non-perturbative parameter ao (/■*/) from % 2 fits of this expression to the mean values of 
thrust, heavy jet mass, C-parameter, total and wide jet broadening. The fits included 
the measured mean values published by various experiments at different centre-of-mass 
energies @, [l9j between 14 and 183 GeV. For the central fit results of a s , we fixed the 



renormalisation scale factor x M to one and the infrared scale to fj,j — 2 GeV. Theoretical 
systematic uncertainties were assessed by varying from 0.5 to 2, fij from 1 to 3 GeV, 
and M. within the quoted uncertainty. 

Figure ^| shows the data values and the fit curves for the C-parameter, the total and 
the wide jet broadening, Bt and Bw- The x 2 /d.o.f. of the fits is between 0.7 (for the 
heavy jet mass) and 1.3 (for the thrust). The numeric results of all fits are summarized in 
Table ^, presenting the values for a s (M z o) and for a (/•£/), their experimental errors and 
systematic uncertainties. We also quote averages of the individual a s and d results which 
were calculated according to the procedure used in Section |1 We consider these results as 
a test of the new theoretical prediction || [7] for the power corrections to the mean values 
of event shapes. It should be noted that the individual results for a (/^/) scatter around 
the average value of a (2 GeV) = 0.473^o;o4i, where the error is the total uncertainty. The 
r.m.s. of the scatter is 0.076. The theoretically expected universality of a® is, therefore, 
observed at the level of better than 20%. 

Our final value of the strong coupling from the fits to the energy dependence of the 
mean values of event shapes is 

a s (M z o) = 0.1177 t° Q S , 

where the error is the total uncertainty. This result lies about 2% above the value obtained 
in our previous analysis || , which did not consider the new Milan correction nor the revis- 
ited calculation for the power correction of the jet broadening observables. Our updated 
result is in good agreement with the world average value ^(J of a™ a '(Mzo) = 0.119±0.004 
and is of similar precision. 



6 Summary 

Inclusion of JADE data and the availability of improved calculations substantially in- 
creased the significance and consistency of QCD tests based on measurements of hadronic 
event shapes. 
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We updated our determinations J| of the strong coupling constant a s at ^/s = 35 and 
44 GeV using ln(i?)-matching for the combination of 0{a 2 s ) and resummed calculations 
now including the C-parameter. Improved perturbative calculations to the jet broadening 
measures Bt and Bw were applied. We found these calculations to describe the data 
better. The a s values were also found to be more consistent with those from other event 
shape observables. Evolving our combined a s measurements to ^/s = M z o we obtain 
0.123 toxtot and 0.123 _oxfol from 35 and 44 GeV data, respectively, resulting in a combined 
value of a s {M z o) = 0.123 ±aool- 

We further investigated mean values of the event shape distributions for the energy 
range of the PETRA and LEP colliders and compared directly with analytic QCD pre- 
dictions plus power corrections for hadronisation effects, the latter involving a universal 
non-perturbative parameter a . Our present studies, based on two- loop calculations for 
the event shapes, yield a s (M z o) = 0.1177 taoo34 5 which is in good agreement with our 
results from the (9(a^)+NLLA fits and also with the world average value. The expected 
universality of the non-perturbative parameter a is now observed at a level of better than 
20%. 
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Tables 



mean value 



1/a ■ da/dC 













35 GeV 












44 GeV 







00- 


-0 


08 





064 


± 





004 


± 





034 





093 


± 





008 


± 





031 





08- 


-0 


12 





434 


± 





019 


± 





079 





801 


± 





047 


± 





131 





12- 


-0 


16 


1 


383 


± 





039 


± 





068 


1 


953 


± 





086 


± 





250 





16- 


-0 


20 


2 


436 


1 

± 





056 


± 





068 


3 


016 


l 

± 





116 


± 





189 





20- 


-0 


24 


2 


678 


± 





060 


± 





229 


3 


172 


± 





123 


± 





288 





24- 


-0 


28 


2 


845 


± 





062 


± 





260 


2 


759 


± 





116 


± 





141 





28- 


-0 


32 


2 


289 


± 





054 


± 





208 


2 


457 


± 





108 


± 





286 





32- 


-0 


36 


2 


O O T 

327 


1 

± 





056 


± 





c\ r* c\ 

269 


1 


748 


1 

± 





086 


± 





098 





36- 


-0 


40 


1 


o o T 

837 


1 

± 





047 


± 





109 


1 


r* c\ o 

623 


1 

± 





082 


± 





o o r\ 

230 





40- 


-0 


44 


1 


441 


± 





041 


± 





131 


1 


163 


± 





068 


± 





150 





44- 


-0 


48 


1 


254 


± 





038 


± 





066 


1 


159 


± 





070 


± 





210 





48- 


-0 


52 


1 


024 


± 





034 


± 





071 





847 


± 





056 


± 





071 





52- 


-0 


58 





839 


± 





025 


± 





034 





720 


± 





041 


± 





070 





58- 


-0 


64 





702 


± 





022 


± 





095 





626 


± 





038 


± 





046 





64- 


-0 


72 





532 


± 





017 


± 





030 





503 


± 





030 


± 





050 





72- 


-0 


82 





453 


± 





014 


± 





051 





348 


± 





022 


± 





034 





82- 


-1 


00 





090 


± 





004 


± 





010 





079 


± 





008 


± 





015 



0.3673 ± 0.0013 ± 0.0040 



0.3404 ± 0.0023 ± 0.0037 



Table 1: Event shape data at s/s = 35 (left) and 44 GeV (right) for the C-parameter ob- 
servable. The values were corrected for detector and for initial state radiation effects. The 
first error denotes the statistical and the second the experimental systematic uncertainty. 





Observable 




stat. 


exp. 


MC 


renorm. 


35 GeV 


Bt 


0.1489 


±0.0008 


±0.0014 


±0.0099 
-0.0048 


±0.0136 
-0.0107 


B w 


0.1367 


±0.0009 


±0.0026 


±0.0089 
-0.0045 


±0.0096 
-0.0075 


C 


0.1480 


±0.0009 


±0.0014 


±0.0097 
-0.0058 


±0.0138 
-0.0110 


44 GeV 


Bx 


0.1458 


±0.0014 


±0.0037 


±0.0072 
-0.0032 


±0.0127 
-0.0100 


B w 


0.1318 


±0.0016 


±0.0051 


±0.0054 
-0.0026 


±0.0081 
-0.0061 


C 


0.1470 


±0.0017 


±0.0027 


±0.0073 
-0.0044 


±0.0133 
-0.0107 



Table 2: Values of a s obtained from fits of 0(a 2 s ) +NLLA QCD calculations to the 
distributions of total and wide jet broadening, B? and By/, and C-parameter, at y/s = 35 
and 44 GeV. Additionally, the statistical (stat.) and the experimental errors (exp.) of the 
fit, the uncertainties due to the Monte Carlo modelling (MC) of the hadronisation and 
due to the choice of the renormalisation scale (renorm.) are given. 
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(a) 


(1 -T) 


(M 2 H /s) 


1 Ti \ 

(B T ) 


/ 7~> \ 

(B w ) 


(C) 


average 


a s (M z o) 


0.1198 


0.1141 


0.1183 


0.1190 


0.1176 


0.1177 


Q range [GeV] 


13-183 


14-183 


35-183 


35-183 


35-183 




X 2 /d.o.f. 


52.2/39 


22.0/33 


22.1/25 


18.8/26 


18.8/16 




experimental 


±0.0013 


±0.0010 


±0.0016 


±0.0020 


±0.0013 


±0.0016 


= 0.5 


-0.0049 


-0.0026 


-0.0038 


+0.0017 


-0.0043 


-0.0027 


x M = 2.0 


+0.0061 


+0.0037 


+0.0048 


+0.0003 


+0.0053 


+0.0026 


M - 20% 


+0.0011 


+0.0013 


+0.0008 


+0.0005 


+0.0009 


+0.0008 


M + 20% 


-0.0011 


-0.0001 


-0.0007 


-0.0005 


-0.0009 


-0.0005 


m = 1 GeV 


+0.0025 


+0.0013 


+0.0017 


+0.0011 


+0.0020 


+0.0014 


\ii = 3 GeV 


-0.0019 


-0.0011 


-0.0014 


-0.0009 


-0.0016 


-0.0012 


Total error 


+0.0068 
-0.0055 


+0.0043 
-0.0030 


+0.0054 
-0.0044 


+0.0029 
-0.0022 


+0.0058 
-0.0049 


+0.0035 
-0.0034 



(b) 


(1-T) 


(Ml/s) 


(Br) 


(B w ) 


(C) 


average 


a (2 GeV) 


0.509 


0.614 


0.442 


0.392 


0.451 


0.473 


experimental 


±0.012 


±0.018 


±0.015 


±0.028 


±0.010 


±0.014 


x M = 0.5 


+0.003 


+0.011 


+0.020 


+0.109 


+0.005 


+0.018 


x M = 2.0 


-0.002 


-0.005 


-0.014 


-0.042 


-0.003 


-0.009 


M - 20% 


+0.058 


+0.084 


+0.046 


+0.032 


+0.050 


+0.053 


M + 20% 


-0.040 


-0.064 


-0.031 


-0.022 


-0.034 


-0.037 


Total error 


+0.059 
-0.042 


+0.087 
-0.067 


+0.052 
-0.037 


+0.117 
-0.055 


+0.051 
-0.036 


+0.058 
-0.041 



Table 3: Values of a s (M z o) (a) and ao(/-tj) (b) derived using \ij = 2 GeV and x M = 1 and 
the 0(a 2 s ) calculations plus two-loop power corrections, which include the Milan factor 
0, and the revisited power corrections for the broadening observables M. In addition, 
the statistical and systematic uncertainties are given. Signed values indicate the direc- 
tion in which a s (Mzo) and ao(fii) changed with respect to the standard analysis. The 
renormalisation and infrared scale uncertainties are treated as an asymmetric uncertainty 
on a s (M z o) and «()(/•*/)■ No error contribution from the infrared scale \xi is assigned to 
a (/i/). 
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Figures 




Figure 1: Distributions for C-parameter (bottom) measured at = 35 and 44 GeV 
and corrected to parton level. The fits of the new 0(a;?)+NLLA QCD calculation for the 
C-parameter are overlaid and the fit ranges are indicated by the arrows. The error bars 
represent statistical errors only. Also shown are the hadronisation correction factors (top). 
The shaded bands represent the uncertainties due to the modelling of the parton shower 
and string fragmention by JETSET 6.3. 
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Figure 2: Values of a s (35 GeV) and a s (44 GeV) derived from C(a^)+NLLA fits to event 
shape distributions. The experimental and statistical uncertainties are represented by the 
solid error bars. The dashed error bars show the total error including hadronisation and 
higher order effects. The shaded region shows the one standard deviation region around 
the weighted average (see text). 
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Figure 3: Energy dependence of the mean values of the total {{B T )) and wide jet broad- 
ening {{By/)), and of the C-parameter ((C)) are shown [19]. The solid curves are the result 
of the fit using perturbative calculations plus two-loop power corrections which include 
the Milan factor |J and the revisited power corrections to jet broadening observablesj7|. 
The dashed line is the perturbative prediction using the fitted value of a s {M z o). 
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